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Abstract 

22 . Solutions of the one dimensional Dirac equation with piece- wise constant potentials are presented 

Z^ \ following standard methods. These solutions show that the Klein Paradox is non-existent and 

represents a failure to correctly match solutions across a step potential. Consequences of this 
exact solution are studied for the step potential and a square barrier. Characteristics of massless 
Dirac states and the momentum linear band energies for Graphene are shown to have quite different 



i-G ' current and momentum properties. 
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I. INTRODUCTION 

The renewed interest in graphene[l| and the close analogy of its band structure to the 
spectrum of the zero mass Dirac equation suggests that a re-examination of several aspects 
of the one dimensional Dirac equation should be carried out. The first of these aspects is 
the well known Klein paradox[2] which continues to persist in the literature. The second 
of these aspects is the question of how closely the graphene spectrum resembles the Dirac 
spectral properties and states. Dragoman in an excellent paper [3|] has recently examined 
some of these issues and has noted that there is no Klein Paradox. In this paper we 
will examine the mathematics of solutions to the one dimensional Dirac equation in the 
presence of a piece- wise constant potential step (the Klein problem). The solutions of this 
differential equation are exact and can be carried out analytically. The usual method of 
analysis will involve the delta function normalization of wave functions in the continous 
spectrum. In order to clarify some of the mathematical details a slightly more meticulous 
method of calculation originally suggested by Von Neumann [J] will be followed so that all 
of the mathematical inferences can be carried out transparently. In the first section details 
of the solutions of the differential equation will be discussed. At the end of that section 
the wave function solutions for the potential step will be exhibited. This will show that 
the Klein paradox arises because of a mis-application of the solutions of piece- wise constant 
potentials for these one dimensiional Dirac equations. In the next section will be a short 
discussion of the admixture of free negative energy states into the positive energy states 
for this step potential system if the potential were to be turned on instantaneously. This 
admixture should monitor the creation of electron-positron pairs near the potential step. 
After this admixture discussion a characteristic of orthogonality in the thermodynamic limit 
will be discussed. The reflection coefficient for finite width square potential barriers will be 
displayed and finally the effect of piece-wise constant potential steps on zero mass states will 
be discussed as well as a comparison with both the energy spectrum and current densities 
for a "one dimensional" Graphene. 



II. SOLUTIONS OF PIECE-WISE CONSTANT POTENTIAL DIFFERENTIAL 
EQUATIONS 

One dimensional solutions of the Sclirodinger and Dirac equations essentially reduce 
to finding the finite, continous, and differentiable solutions of a Sturm-Liouville type of 
differential equation on the whole line interval (— C)0, oo)|5|,[7|. For the discrete spectrum 
of such equations there is an orthonormal set of eigenstates. For the continuous spectrum 
the momemtum eigenstates are functions that can only be normalized with Dirac delta 
functions. The delta function normalization is generally quite satisfactory except for some 
calculations involving squares of delta functions. 

Because of a desire to be mathematically careful and transparent in this examination of 
the Klein paradox, a procedure will here be followed that dates back to von Neumann[4] 
and many others[6i]. Specifically, all of the wave functions to be used here will be "box 
normalized" in regions of length L and final results will be obtained in the thermodynamic 
limit {L -^ oo). Within the interval of length L the allowed momentum values will be 
discrete k = 2TTn/L, where n is an integer. This whole process will not be essential for 
the demonstration that the Klein paradox is not a paradox at all, but will be useful in 
studying the admixture of free particle states and those interacting states in the presence 
of the potential step. For each finite L we will find that the eigenvalue spectrum will be 
discrete and that in the thermodynamic limit (L -^ oo) the spectrum becomes continuous 
and the eigenstates become delta function normalized. Details of the transition to the 
thermodynamic limit have been discussed in detail by Messiah 



Sneddon 



71], Arfken and Weber [8|], and 



9| and some of these details will be discussed in the following. 
The Dirac equation with a constant potential has exact solutions which are the same as 
the free particle solutions except that the energy Ek can be different from the free particle 
case 



Ek = Vo± Vm2 + k^ (1) 

by the addition of the constant potential. In order to simplify the notation in this study 
the velocity of light c = h = 1 will be used in the formulas to follow. The Dirac equation 
will be studied with 2x1 spinors for simplicity. This is a standard simplification sometimes 
described as "no spin" or "single spin". The basic idea of this paper is that we consider 
these solutions to Dirac equations in two regions {—L, 0) and (0, L) and construct all wave 



functions out of two parts 

^{z)=MzM-z)+^r{z)Q{z) (2) 

where Q{z) is the Heaviside function (unit step function at zero), iPl{z) is a solution in the 
left side interval (— L, 0) and iPr{z) is a solution in the right interval (0, L). The boundary 
condition is that the two solutions be continuous at the origin 

M^~) = M^^)- (3) 

In each region, (— L,0) and (0,L),the orthonormal positive and negative energy states, 
respectively, are given by 




and 

where 

Uk = . (6) 

m + \/m'^ + A;2 

Both of the orthonormal solutions above are written for a positive-direction probability 
current density Jk which is given by 



2u 



k 



Reversing the k value gives the current in the negative direction carrying states. The 
energies of these states will depend on the side of the origin. On the left the energy will be 
±vm^~+~P and on the right Vq =t \/m? + /c'^. In the following paragraphs the wave vectors 
on the right side {z > 0) will be indicated with a prime as k'. 
On either side of the origin these solutions are orthonormal 

pi{k—k')L 1 

{nk\nk') = {pk\Pk') = ^k,k' = .., _ ,,.j- (8) 

where Ak^k' is the Kronecker delta, Ak^k' = 1 if k—k' and zero if {k — k')L = 27r{n—n') ^ 0. In 
each interval the allowed momentum values are k = 27rn/L, where n is any integer. Similarly, 
{nk\Pk') = 0. The discrete k values are those which guarantee that the momentum operator 



is hermitian, namely, that ■?/''^(0)-?/'(0) = ip^ {L)ip[L) for the right side, and similarly on the 
left side. The usual delta function nromalization is achieved by removing the normalization 
factor 1/vZ from each wave function and recalling that for all sequences leading to the 
thermodynamic limit 

lim LAfc,fc' = 5{k - k') (9) 

where 5{k — k') is the Dirac delta function. In the following it will often be the case that the 
normalizations for each half interval 1/(vLa/1 + m^) will be discarded and normalization 
on the whole interval {—L^ L) will be re-calculated. 

For the Dirac equation on the whole interval (—00, 00) there are no states in the mass gap 
—m < E < m. But, for each of the two half intervals there are evanescent (exponential) 
wave functions that preserve momentum hermiticity. These are labeled by a complex 
momentum in and have both positive and negative energy relative to the center of the mass 



gap Vo ± y/m'^ — K^ and can be written for each fixed L in the sequence as 



0^>+)(z)=\ ^ )e^''" + e-"^[ ^ I e''" (10) 



Oi>~\z) =1 '"^M e-- + e-'^^ '""^ I e- (11) 



1 



where 



y^^ = ^ A — r (12) 

m + v'^ — t^ 
for the interval (0, L). Similar functions Oi and Oi are defined for the interval {—L, 0). 
In the thermodynamic limit the second term vanishes and these functions yield the usual 
single exponential function. 

When the step potential Vq is zero, there are no allowed states in the mass gap. Even 
when < Vo < 2m there are no states in the part of the gap which remains on both sides 
of the origin. 

The general solution is now achieved by equating at the origin the two functions, one 
from from each side at the origin, that have the same energy, and using this equality of 
energy to determine the relationship between the momenta on the two sides. In Tables (1 
and 2) are shown the energy ranges and the structure of the wave functions in each region 
which have the same energy. Note that for simplicity, we are examining states that initiate 



with a current from the left and contain a reflected wave on the left and transmitted wave 
on the right. 

There are two distinct Vq energy ranges and to simplify the discussion Table 1 shows the 
wave function pairings and the relationships between the left and right side momenta for 
small < Vq < 2m and Table 2 shows the same for Vq > 2m. The table shows the wave 
functions and the relationship between the momenta on each side of the origin and for each 
energy range. 

Referring to Table I, for small Vq (0 < Vq < 2m,), we will discuss the states from the 
top of the Table (highest energy) to the bottom (lowest energy). At the top are the 
states for which E^ > Vq + m. As can be seen in the Table we have "positive energy" 
functions Pk{z) + fk P-k{z) on the left (/^ is the reflection amplitude) and on the right we 
have QkPk'iz) where gk is the transmission amplitude. The next set of functions are in the 
energy interval m, < Ek < Vq + m,, which are of the form Pk{z) + e^'^p-k{z) in the left side 
and are matched with an evanescent wave 0^{z) on the right side. Because the evanescent 
states carry no current, the reflection amplitude has unit magnitude and thus gives complete 
reflection. The next energy range Vq — m < E < m corresponds to the empty mass gap 
with no wave functions satisfying the boundary condition. The next lower energy range 
—m < Ek' < Vo — m, corresponds to another combination of which half is evanescent and 
cannot carry any current. These states have an evanescent wave 0-K.{z)on the left and 
nk{z) + e^'^n_k{z) on the right. The reflection amplitude again has magnitude 1 as found 
before for a combination of traveling and evanescent waves. The lowest energy states are 
the fully negative energy states and have n^k{z) + fk nk{z) on the left and gkU^k'^z) on the 
right. 

As the potential step size Vq increases and enters into the range Vq > 2m,, a new pairing of 
functions that Klein did not consider begins to appear. These states are discussed Table 2. 
Examining the top of Table 2, the highest energy range Ek > Vq + m is unchanged from the 
previous discussion. The second energy range is now expanded to Vq — m, < E,^ < Vq + m, 
and remains, as before, with complete reflection. The new range m < Ek < Vq — m 
pairs wave functions Pk{z) + fk P-k{z) on the left with transmitted functions gkf^-k'{z) on 
the right. The remaining lower energy ranges and functions are not changed significantly 
from the previous discussion. This new set of solutions on both sides in the energy range 
mi < Ek < Vq — m is what Klein did not consider. In the following we will solve for 
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the reflection and transmission amplitudes for those cases with non-zero probability current 
density. 

For energies above the step, Ei^ > Vq + m it is straightforward to show that 

Vm2 + A;2 = Vo + \/m^ + A;'2 (13) 



(Wfc - M^) _ 2m, 

Jfc — 

where 



•^*^ ~ 7 ^ — T\ ' 5''= ~ 7 \ — M ' (1^) 






and 



, , ^w? + k'^ -Vo-m 

""^ = V / 2^^2 T/ "^- ^16) 

y/m^ + k^ — Vq + m 

Using the usual definitions of current density, it is easy to determine that the reflection 

coefficient R and the transmission coefficient T are given by 

and 

from which it is easy to see that R + T = 1. 

For the ranges of energies in which the waves on the right (or the left) are evanescent, 
we easily find that T = and R = 1. 

For the case which Klein did not consider, when Vq > 2m similar calculations yield 

Vm2 + k^ = Vo- Vm^ + A;'2 (19) 

_ {uku'f, - 1) _ 2uk 



^k=\l / 2^,2^ • (21) 

Vq — yvn^ + k"^ + m 



and Uk is defined as before. The refiection and transmission coefficients are: 



R=\h\' = P^^^ (22) 



and 



_ < 2 _ ^UkU'^, 



from which it is also easy to see that R + T = 1 . 

From these calculations it is clear that there is no paradox. Klein simply did not match 
the appropriate solutions in the two regions. This mis-match is what is responsible for the 
usual assumption of a failure of particle conservation (usually taken to imply the production 
of particle hole pairs near the potential step) that has been so often interpreted as the 
meaning of the Klein paradox. These analytic solutions show that particle number is always 
conserved independent of the size of the step potential. The surprising phenomena is that 
the barrier is close to being transparent in a certain energy range if the step height Vq is large 
enough. This behavior of the reflection coefficient is shown in Fig. 1 where the reflection 
coefficient R has been simultaneously plotted for three different values of Vq = 0.5m, 3m, 
8?72. In all of these cases there is a region of complete reflection with an energy width of 2m 
and centered on the value of Vq. At higher energies the reflection coefficient decreases with 
increasing energy. 

The surprising result is that the reflection coefficient decreases in the energy interval 
fn ^ Ek < Vq — m. and the barrier becomes partially transparent in this energy range. This 
result is an exact consequence of the matching of exact solutions in the two regions joined by 
the boundary condition. On physical grounds this result would seem to be unexpected, but 
the potential step has pulled what used to be negative energy states into a positive energy 
range and their ability to carry current leads to the partial transparency of the potential step 
in this energy range. One possible origin for this surprising phenomena could be with Dirac's 
original choice for adding potential energies onto the free particle Dirac Hamiltonian. 

III. OVERLAP BETWEEN STATES OF NON-INTERACTING AND INTERACT- 
ING HAMILTONIANS 

Another aspect of the behavior of this system can be probed by studying the overlap of 
these interacting states in the presence of the potential step with the free states of the free 
particle Hamiltonian. This is equivalent to asking how the wave functions are matched in 
the sudden approximation if the step potential were instantaneously turned on. If there 
were a production of extra electron positron pairs by the potential step, it could be expected 
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that this overlap would be an indicator of such when the potential Vq becomes larger than 
2m. In order to carry out this estimate we would need to calculate 

N=J2 IK" in) I', (24) 

k",Pk 

where nk"{z) is a negative energy wave function of the zero potential Hamiltonian spanning 
both intervals (— L,0) and (0,L), and Pk{z) represents any positive energy state of the 
Hamiltonian with the potential step Vq. And the expression is summed over all positive 
energy states with non-zero matrix elements. 

It is first important to notice that all of the positive energy states of on both sides of the 
potential step, which are above the mass gaps on each side, will have no contribution to this 
matrix element because all of these states are made up of linear combinations oi pk{z) and 
p^k{z) which are orthogonal to nk"{z) in both of the intervals (— L, 0) and (0,L). So, the 
only states that can overlap with the negative energy states can be: 

(1) evanescent states in the mass gaps (which will be shown to be negligable in the 
thermodynamic limit, see Appendix I), 

(2) the negative energy states nk{z) + e^'^n-k{z) in the range —m < Ek' < Vq — m which 
are matched with evanescent states. The wave functions for these are of the form 

^2^ = A2{{0^^{z)Q{-z) + (rikiz) + e''^n.k'{z))Q{z)), (25) 

and 

(3) positive energy states appearing when Vq > 2m which carry current on both sides 

^3(^) = A^iipkiz) + fkp-k{z))e{-z) + {gkn.k'{z))e{z)). (26) 

Before evaluating these three cases, it is possible to arrive at an intuitive estimate of 
this quantity by simply asking from a density of states perspective how many of these 
originally negative energy states have been pulled up from energy — -m to positive energies 



Vo — ^/mF+W by the magnitude of Vq. If we ignore the boundary conditions, the number 
of such states would be given by the integral 



L /'v'^»^^"+^")^,_ Lv/\/o(\/o + 2m) 



27r Jo 27r 

Let us first examine case (2). The first step is to evaluate the normalization A2 which 

yields 

A,{k) = ^ ^ ^==. (28) 



Neglecting terms of order 1/L and smaller and using the orthogonality of nk'{z) and nk"{z) 
it is straightforward to show that 

{nk"\-^2ik)) = -(Afc»,fc, + e''f'Ak",-k'). (29) 

Now, for each L in the sequence to the thermodynamic limit, it is true that 

^k",k' = ^k",k' (30) 

so we obtain 

|(nfc"|^2(A;))|' = ^(Afc",fc. + Ak",-k'). (31) 

This implies that 

N2 = 5^ |(nfe.|^2(A;))P = ^ 5^ 1 = ^VVo{Vo + 2m) (32) 

k",k k' 

for the energy range —m < Ej./ <Vq — m. If Vq > 2m, the lower limit of the k' integration 



becomes v^Vo(Vo — 2m) so that the expression for N2 in that case becomes 

N2 = ^( VK)(K) + 2m) - y/Vo{Vo - 2m)). (33) 

4:71 

Now let us examine case (3). In this case the matrix element becomes 

{nkii\^^{k)) = A'i{k)gk I nk"{z)^n_k'{z)dz (34) 

Jo 



where 



Asik) = , (35) 



Substituting for /^ and gk from the eqn. (20) and (21) above, we find 



^ulul,{3 + ul) + {1 + 3ul) 



So, the contribution from this case is 



L fV^o(Vo-2m) 2ulil+ul,] 



iVg = — / ^^^^- ^•\; f^ :^dk. (37) 

2nJ, ulul,{3 + ul) + {1 + 3ul) ^ > 

The remaining case 1 is analyzed in Appendix I and gives a result which is not extensive 
with the length L and so makes no contribution in the thermodynamic limit. The details 
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and certain aspects of orthogonality in the thermodynamic hmit are presented in Appendix 
I. 

The dependence of A^ on the barrier height Vq is shown in Fig. 2. For small Vq the 
intuitive estimate and the exact result agree. When the potential height becomes greater 
then 2m, the lower limit of Eqn. (33) and the integral for N^ makes their contribution. The 
smooth curve is the intuitive estimate given by Eqn. (32) and the lower curve with the kink 
at Vq = 2m is the exact result for this quantity. Note that the behavior of this quantity is 
smooth after the threshold Vq = 2m and less than the intuitive estimate. 

IV. TRANSMISSION THROUGH SQUARE BARRIERS 



Walker and Gathright lO] worked out all possible one dimensional transfer matrices for 
the non-relativistic Schrodinger equation across potential discontinuities. They constructed 
transfer matrices for any arrangement of potential discontinuities by building the transfer 
matrices out of products of two different matrices: the discontinuity matrix d and the 
propagation matrix P. All such possible matrices have been worked out for the Dirac 
equation and will be presented elsewhere. For the one dimensional Dirac equation the 
discontinuity matrix d is of the form 



d{a, h) 



a + b a — b 
a — b a + b 



(38) 



and the parameters a and b are replaced by 1 or by u^ or iw^ or by ratios of these depending 
on the energy range in the segment between two discontinuities of the potential. The 
propagation matrix is identical to those found by Walker and Gathright 



P(a) 



e" 
e"" 



(39) 



Once relativistic matrices d and P have been constructed for all possible pairs of states 
which are possible at a discontinuity, it is immediate to construct the transfer matrices that 
will connect the two solutions at the ends of a region / where the potential has a constant 
value Vi. A simple, symmetric example of the use of these transfer matrices is a square 
barrier whose height is Vq = 5.5 and whose width is a/L = 5/150 = 1/30. A plot of the 
reflection coefficient R versus the momentum of the initial wave k is shown in Fig. 3 . 
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Notice the Rainsauer ininima (peaks in T ) in the Reflection coefficient. Ill] Notice also 
that the Reflection coeflicient for this symmetric case behaves much hke the step potential 
if the potential energy Vq is larger than the mass gap width and becomes semi-transparent 
at lower energies. 



V. PROPERTIES OF ZERO MASS SOLUTIONS OF THE DIRAC EQUATION 



The well known analogy between the band structure of graphene and the energy spectrum 
of the massless Dirac equation has received much attention recently. For the non-relativistic 
band structure calculations the velocity of a band state represented by the energy e^ is linear 
in \k\ and is given by e^ = ±f |/i;| (in one dimension). For electronic band structures the 
current density is proportional to the group velocity which is given by the derivative of the 
band energies, 

and the :" graphene" band structure spectrum will contain both positive and negative veloc- 
ities for both positive and negative energies (relative to the center of the band). 

Examining the zero mass eigenstates of the Dirac equation shows that there are some quite 
different behaviors than those observed in the non-relativistic band structure for Graphene. 
The (1+1) Dirac equation has the following simple form in a constant potential Vq 



^ 



idz 




Akz 



{E-Vo 




Akz 



(41) 



has eigenstates 



^±{k) 



V2 



1 

±1 



ikz 



with correspond to the energies 



E-Vo = ±k. 
The probability current densities carried by these states are: 

1 



J4 



^l(T,;*_ 



and 



J_ 



1 1 



\&^ a,.\&_ 




-1. 



(42) 



(43) 



(44) 



(45) 
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So the positive direction currents are carried only by positive energies E — Vq > and 
negative direction currents are carried hj E — Vq < 0. The massless Dirac equation does 
not allow negative currents to be carried by positive energy eigenstates. 

If we consider the reflection at a potential step aX z = and apply the boundary condition 
we obtain, 

k = Vo + k' (46) 



and 





9k' \ . (47) 



which yields the equations 

l + fk = 9k' (48) 

l-fk= 9k' (49) 

which implies that fk = ^ and Qk' = 1. This condition implies that the potential step 
makes no reflections for a zero mass particle. This condition and the fact that current 
direction is so strongly associated with the sign of the energy (relative to Vq ) indicates that 
the analogy between the graphene bandstructure and the massless Dirac equation is not 
completely accurate. 

VI. CONCLUSIONS 

The Klein Paradox is not a paradox. It is simply a mis-application of the processes by 
which the solution of piece-wise constant potential differential equations are constructed. 
When the appropriate wave functions at the same energies are connected, the reflection and 
transmission coefficients are continuous functions of the incident wave vector and always 
obey the conservation of particle number. The surprising characteristic of these solutions 
is the near transparency of the step potential at low energies if Vq > 2m. This property of 
the solutions arises from the fact that for large Vq states which were originally at negative 
energies are now pulled up into positive energies and it becomes possible for a current to 
be carried through the step. This result is clearly a property of the solutions of this Dirac 
equation. The more difficult question is whether this behavior is physically to be expected. 
This property reflects a choice made by Dirac when he decided to add a potential energy 
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to his free particle equation. He chose to add it to the ex. ■ p as opposed to adding the 
potential ener gy to the mass m. These questions have been examined to some extent in 



other contexts 



121. 



The connection between the overlap of the negative energy states of the non-interacting 
Hamiltonian and the Hamiltonian with the potential step, if the latter is turned on instan- 
taneously was examined by studying the summation of the overlap matrix element between 
positive energy states and the initially negative energy states. It was shown that the 
accurately determined overlap was consistent with an intuitive picture that the overlap 
represented the number of negative energy states that have been pulled above the energy 
E = —m. As should be expected because the Klein paradox does not exist, there is no 
anomalous behavior of this overlap as the potential step is increased above the threshold 
Vo = 2m. 

Generalizing the transfer matrices from the non-relativisitic Schrodinger equation to the 
Dirac equation allows the treatment of a variety of potential barriers and steps and, numer- 
ically, any smooth potential that can be approximated by piece-wise constant potentials in 
short intervals. By way of an example, the case of a square potential barrier was briefly 
discussed. The presence of Ramsauer resonances and the transparency of the barrier were 
found in direct analogy to the results of the step potential. 

Finally, zero mass eigenstates of the Dirac equation were examined. Positive current 
densities were only carried by positive energy states (relative to Vq) and negative current 
densities were only carried by negative energy states. It appears that a step potential at 
the origin has no effect on these states. Both of these conditions are quite different from 
the band structure of Graphene which has stimulated the analogy between that material 
and the solutions to the Dirac equation. The failure of a step potential to influence the 
zero energy states seems to be quite unphysical, and seems again to be related to Dirac's 
original choice by which he added the potential energy to the free particle Dirac equation. 
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VIII. APPENDIX I: OVERLAP BETWEEN EVANESCENT AND FREE NEGA- 
TIVE ENERGY STATES 

In this appendix the overlap of the negative energy free particle eigenstates with the 
evanescent states in the gaps are examined and found to be of order 1/L. 

As an example, consider a state where the overall energy is in the range Vo < E < Vo + m. 
One of the wave functions in this range is 



where 



^,{z) = A,{k){{pk{z)+e%.k{z)M-z) + {g.e-^i . )e(^)), (A-1) 



A,{k) = , ^ ^^ . (A-2) 

The matrix element, ignoring terms of order e~'^^, is 

{nk"\'^i{k)) = ^ • (A-3) 

^ v^v/(l + OV2(l + ul) + ^il + wl){K - ik") 

The contribution of this matrix element to A'^i as L ^ oo is given by 

-. - g l(-l*.)P . ^/*"/..^^^-j|i±|,,^. ,A.4, 

Since this integral is not extensive in L when compared to the other contributions which are 
proportional to L, this contribution is negligible in the thermodynamic limit. 

The fact that the overlap olnkii{z) with the evanescent states Oi^{z) was negligable in the 
thermodynamic limit raises a question involving evanescent states in the mass gap and more 
generally parts of wave functions within the confines of a square potential barrier Q < z < 
a < L. This detail has already been discussed by Arfken and Sneddon in their discussion of 
the evolution of the fourier integral from the Fourier series in the thermodynamic limit. In 
their discussion it is observed that in the transition L — > oo the constant term ao/2 of the 
Fourier series becomes negligible in the thermodynamic limit. 

Similarly, the corresponding relationship between "box normalized" wave functions in 
the thermodynamic limit becomes apparent in the theorem for Hermitian operators that 
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eigenvectors with different eigenvalues must be orthogonal. Consider two states of the type 
from case 1 with different wave vectors, and thus different energies 

^u,{z) = A,{h){{p,M + e''^p^,,iz))Qi-z) + ig.,e-^^'(. ^ ))Q{z)) (A-5) 



^,,{z) = A2{h){{pkM + e'^P-k,{z)M'^) + {9..e~^'^. ^ ))e(^)) (A-6) 

and the energies are: 



Eu, = \Jm^ + kl ^ s^m^ + hi = Ek,. (A-7) 

If we evaluate the overlap between these two vectors, the integrals on the left hand interval 
(— L, 0) immediately give zero because of the orthogonality of the pk^ functions. This leaves 
an integral on the interval (0, L) 

{^kA^k.) = A,A,iO+ ,^pP'^l^-^'''\ J - (A-8) 

On first refiection this second term is not zero and appears to violate the orthogonality 
theorem, but it must be noted that the factors A1A2 provide a factor of 1/L and this matrix 
element vanishes in the thermodynamic limit. The orthogonality of the kind of states 
combined in the process of joining solutions of piece-wise constant potentials is dominated 
by the momentum eigenstates and integrals over smaller intervals can be non-zero for finite 
L, but make no contribution in the thermodynamic limit. 

IX. REFERENCES 



[1] K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Katsnelson, I. V. Grigorieva, S.F. 

Dubonos, A. A. Firsov, Nature 438, 197(2005), S. Y. Zhou, G.-H. Gweon, J. Graf, A.V. 

Fedorov, C.D. Spataru, R.D. Diehl, Y. Kopelevich, D.-H. Lee, S.G. Louie, A. Lanzara, Nature 

Physics 2, 595-599(2006). 
[2] O. Klein, Z. Phyzik, 53, 157(1929), C. Itzykson and J. Zuber, Quantum Field Theory, McGraw 

Hill, New York, 1980. 
[3] D. Dragoman, ArXiv quant ph/0701083. 



16 



[4] J. Von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton University 
Press, Translated by R.T. Beyer, 1955. 

[5] E. Merzbacher, Quantum Mechanics, Wiley, NY, 1970, Chapt. 6. 

[6] P. Morse and H. Feshbach, Methods of Theoretical Physics, McGraw-Hill, New York, 1963, M. 
Born &: K. Huang, Dynamical Theory of Crystal Lattices, Oxford, London (1954), M. Born &: 
T. von Karman, Phys. Z. 13, 297 (1912). 

[7] A. Messiah, Quantum Mechanics, North Holland, 1961, Vol. 1, Chapter 5 

[8] G. B. Arfken and H. J. Weber, Mathematical Methods for Physicists, Elsevier, 6th ed., 2005, 
p. 936 fF. 

[9] I. N. Sneddon, Fourier Transforms, Dover, New York, 1995, Sect. 3.2. 
[10] J. S. Walker and J. Gathright, Am. J. Phys. 62, 408-422(1994). 

[11] R. L. Liboff, Introductory Quantum Mechanics, Holder, Oakland, CA, 1980, N. Zettilli, Quan- 
tum Mechanics, Wiley, NY, 2001. 
[12] W. Greiner, B. Muller, J. Rafelski, Quantum Electrodynamics of Strong Fields, Springer Ver- 
lag, Berhn, 1985. 

X. TABLES 



Table 1: A tabulation of the energy range and the type of wave functions that are matched 
at the origin as well as the equation relating the wave vector k on the left and k' on the right 
for small Vq < 2m. In this example, the two mass gaps have a significant overlap which 
contains no wave functions. 



E range 


M^) 


M^) 


k and k' 


E > Vo + m 


Pk{z) + fkP-kiz) 


9kPk'{z) 




yjm? + A;2 = I/q + Vm2 + A;'2 


m<E<Vo + m 


Pkiz) + e^'^p^kiz) 


9kO^{z) 




v/m2 + A;2 = Vo + Vm2 - /t2 


Vo — m < E < m 


no w.f. 


no w.f. 


True Mass gap; no states 


-m< E <Vo-m 


9-.0^.{z) 


Ukiz) + e'^^n-k'iz) 




±Vm2 -K^ = Vq- yjm? + A;'2 


E <-m 


n^k{z) + fk nk{z) 


9kn-k'{z) 




-Vm2 + A;2 = Vo - Vm2 + A;'2 



Table 2: A tabulation of the energy range and the type of wave functions that are matched 
at the origin as well as the equation relating the wave vector k on the left and k' on the 
right for large Vq > 2m. The two mass gaps have a significant separation and the function 
pairs that Klein ignored are included in this energy range. 
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E range 


M^) 


M^) 


k and k' 


E >Vo + m 


Pk{z) + fk p-k{z) 


gkPk>{z) 




Vm2 + k^ = Vo + Vm^ + k'^ 


Vo< E <Vo + m 


Pfc(2;) + e^VfeW 


g.O^>t\z) 




\/w? + A;2 = \/o + Vm^ - k^ 


Vo-m<E<Vo 


Pk{z) + e^^p.k{z) 


g.0i-\z) 




y/m"^ + k"^ = Vq — y/m'^ — k^ 


m < E <Vo- m 


Pk{z) + fk P-k{z) 


gkn-k'{z) 




Vm2 + A;2 = \/o - Vm2 + A;'2 


0< E <m 


g-.oi^^z) 


nk'{z) + e''^n_k'{z) 




+Vm2 -K^ = V^- Vm2 + A;'2 


-m< E <0 


g-.oi-^i^) 


nk'iz) + e''f'n.kiz) 




-Vm2 -K^ = Y^- Vm2 + A;'2 


E <-m 


n^kiz) + fk nk{z) 


gkn^k'iz) 




-Vm2 + A;2 = Vo - Vm2 + A;'2 



XI. FIGURE CAPTIONS 

Fig. 1: Three independent plots of the reflection coefficient R versus the wave vector k 
of the incident wave for three different values of Vq = 0.5, 3.0,8.0 rest masses are plotted 
together. Each curve can be identiffed by the value of k at the center of the R=\ plateau. 

Fig. 2: Two curves, one approximate and one exact, for the overlap A^ between the 
positive energy states in the presence of the step potential and free negative energy states 
in the absence of the step potential. The smooth curve is the "intuitive" A^ derived from 
the density of negative energy states shifted up by the potential step. The lower curve 
with the break at Vq = 2 is the complete calculation of A^ as a function of Vq for states 
satisfying the boundary condition. Note that the exact curve does not indicate excessive 
overlap (electron-positron) pairs above the threshold Vq = 2m. 

Fig. 3: A plot of the reffection coefficient R versus the inicident wave momentum k for a 
square barrier with height Vq = 5.5 and width ajL = 1/30. Notice the prominent Ramsauer 
minima. 
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